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Double-polarization asymmetries for inclusive ep scattering were measured at Jefferson Lab using 
2.6 and 4.3 GeV longitudinally polarized electrons incident on a longitudinally polarized NH3 target 
in the CLAS detector. The polarized structure function gi(x,Q 2 ) was extracted throughout the 
nucleon resonance region and into the deep inelastic regime, for Q 2 = 0.15 - 1.64 GeV 2 . The 
contributions to the first moment Fi(Q 2 ) = J gi(x,Q 2 )dx were determined up to Q 2 = 1.2 GeV 2 . 
Using a parametrization for g± in the unmeasured low x regions, the complete first moment was 
estimated over this Q 2 region. A rapid change in Ti is observed for Q 2 < 1 GeV 2 , with a sign change 
near Q 2 = 0.3 GeV 2 , indicating dominant contributions from the resonance region. At Q 2 = 1.2 
GeV 2 our data are below the pQCD evolved scaling value. 

PACS numbers: 13.40.Gp, 13.60.Hb, 14.20.Dh 



Electron scattering has played a long and distinguished 
role in the study of nucleon structure. Inclusive deep 
inelastic scattering (DIS) studies revealed the parton 
constituents (quarks and gluons) of the nucleon. At low 
values of the four momentum transfer Q 2 the nucleon 
structure functions depend on both Q 2 and on the energy 
transfer v = Eq — E' , where Eq and E' are the initial and 
final electron energies. In the asymptotic limit, Q 2 — ► 
00, the scattering is described by perturbative Quantum 
Chromodynamics (pQCD) as the absorption of a virtual 
photon on a single free quark. In this limit the structure 
functions depend only on the Bjorken scaling variable, 
x = Q 2 /2Mv, where M is the nucleon mass. 

When both the incident electron and the target nu- 
cleon are polarized (double-polarization) the scattering 
cross section depends on two additional functions of Q 2 
and v, <?i and 172- In the framework of pQCD these also 
depend only on x as Q 2 — » 00. In this limit gi(x) has a 
simple interpretation as the sum over the x-dependent 
polarizations of the various quark flavors. A double- 
polarization experiment at CERN showed that in the 
framework of asymptotic QCD only a fraction of the nu- 
cleon spin could be attributed to the intrinsic spin of 
the quarks. These measurements were confirmed by sub- 
sequent experiments at CERN, SLAC, and DESY, (see 
the reviews |2j, |3| and references therein). Corrections 
(higher-twist terms) for finite Q 2 result in excellent fits 
to the extensive data set for Q 2 down to the order of 2 
GeV 2 . It is generally accepted that the intrinsic spin of 



the quarks accounts for about 25% of the nucleon spin, so 
that other degrees of freedom, such as gluons and quark 
orbital angular momentum, must account for the rest. 

Until recently, only a few double-polarization exper- 
iments were carried out with energies below 25 GeV, 
an important example being the ground-breaking ex- 
periment at SLAC in the late 1970's. More re- 
cently, the spin structure function g\ and its first mo- 
ment Fi(Q 2 ) = J gi(x,Q 2 ) dx have become a focus at 
lower Q 2 and in the resonance region in experiments at 
SLAC |3, Hermes 0, 0, and JLab H to obtain a 
better understanding of QCD in the confinement regime 
and to study the phenomenon of duality between the res- 
onance region and the deep inelastic region. The rapidly 
changing helicity structure of some resonances as a func- 
tion of Q 2 is expected to have a strong influence on g\. 
Double-polarization measurements at lower energies can 
therefore indicate where pQCD breaks down and deter- 
mine where multi-parton processes and coherent effects 
due to nucleon resonance transitions are important. 

For inclusive scattering of electrons and protons po- 
larized along the axis of the electron beam, the double- 
polarization asymmetry is given by: 
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where cr* and are the cross sections for the electron 
and proton spins parallel and anti-parallel, respectively. 
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The factor e = [1 + 2(1 + v 2 /Q 2 ) tan 2 (6»/2)]- 1 is the vir- 
tual photon polarization for electron scattering angle 9. 
The parameter R = <tl/&t is the ratio of the absorp- 
tion cross sections for longitudinal and transverse virtual 
photons. The kinematical factor -q = e^/Q 2 / (Eq — e E'). 
The angle between the virtual photon, 7* , and the beam 
direction is given by # 7 . 

The photon asymmetries A\ and A2 can be written in 
terms of the virtual photon absorption cross sections as: 
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in which ct^/ 2 and cr^/ 2 are transverse cross sections and 
<jt is half their sum. The superscripts denote the helic- 
ity of the j*p system. The cross section glt arises from 
longitudinal-transverse interference. In the nucleon res- 
onance region the relative contributions of helicity 1/2 
and 3/2 vary from one resonance to another and depend 
strongly on Q 2 . The spin structure function g\ is linearly 
related to A\ and A 2 by: 
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where F\ is a structure function appearing in the unpo- 
larized electron scattering cross section. 

If the contribution from elastic scattering (x = 1) is 
excluded, Ti vanishes at Q 2 = 0, where, in addition, its 
slope is constrained to be ne gati ve by the Gerasimov- 
Drell-Hearn (GDH) sum rule Hi EH for absorption of 
real photons. The validity of the GDH sum rule has been 
experimentally tested to better than 10% by experiments 
at Mainz [i^ and at ELSA On the other hand, Ti 
is known to be positive at high Q 2 [2- [3] Therefore, the 
constraints near Q 2 = and at high Q 2 imply that Ti 
must change sign at some low Q 2 , where it is expected 
to be dominated by the nucleon resonances 0, 0] . 

At low Q 2 , meson-baryon dynamics have been treated 
in chiral perturbation theory by Ji and Osborne |l6| and 
by Bernard, Hemmert, and Meissner 17] to evolve Ti to 
non-zero values of Q 2 . Badelek, Kwiecihski, and Ziaja 
|l8| have used the generalized vector dominance model 
and GDH to make predictions for g\{x) at low Q 2 . At 
high Q 2 the asymptotic value of Ti has been evolved 
down to Q 2 — 1 GeV 2 by using pQCD and the Operator 
Product Expansion 0. Phenomenological approaches, 
either with explicit inclusion of resonance parameters 
|20l l2ll] or with general parametrizations of structure 
functions 22], have been used to cover the entire Q 2 
range. 

The present measurements were carried out with 2.6 
and 4.3 GeV longitudinally polarized electrons incident 
on a longitudinally polarized target located at the cen- 
ter of the CLAS detector [23}. A toroidal magnetic 



field, symmetric about the beam axis, is generated by 
six super-conducting coils. The coils separate the detec- 
tor into six independent spectrometers that use wire drift 
chambers for track reconstruction, scintillation counters 
for time-of-flight measurements, threshold gas Cerenkov 
counters, and lead-scintillator electromagnetic calorime- 
ters. Electrons can be detected and identified for mo- 
menta down to 0.35 GeV/c and for polar angles from 
about 8° to 50°. The polarization of the beam, which 
was measured frequently with a M0ller polarimeter, was 
typically 70%. The beam helicity was flipped at a rate of 
1 Hz in a pseudo-random sequence to minimize system- 
atic effects. 

A microwave pumped solid nuclear target 1241 using the 
method of Dynamic Nuclear Polarization |25j was built 
for the CLAS. Its design is similar to those employed in 
experiments at SLAC |5j. Horizontal cylinders, ~ 1 cm 
long, were packed with ammonia pellets (either 15 NH3 
or 15 ND3). The cylinders were mounted with their axes 
along the beam line on a movable ladder, along with an 
empty cup and a 2.2 mm thick carbon disk. The tar- 
get cell was immersed in a liquid He bath at T=1.2 K. 
A pair of super-conducting Hclmholtz coils coaxial with 
the beam generated a uniform 5 T magnetic field in the 
target volume. A maximum proton polarization of 70% 
was obtained. During beam irradiation the proton po- 
larization typically dropped gradually to 40%, at which 
point the target material was removed and replaced in 
order to restore the original polarization. 

Charged particles were identified from the measured 
momentum and time-of-flight. The Cerenkov counters 
and calorimeters were used to reduce the 7r _ contamina- 
tion in the electron sample to less than 1%. The scattered 
electrons were binned according to Q 2 and the invariant 
mass W = (M 2 + 2Mv — Q 2 ) 1 ^ 2 of the recoiling hadronic 
system. For each bin the polarization asymmetry defined 
in eqn. ^is: 
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in which P e and Pt are the averaged polarizations of the 
electrons and protons, respectively; and Am are the 
number of observed electrons normalized to the total in- 
cident beam flux for parallel and anti-parallel beam and 
target spins; and N^kg is the number of electrons scat- 
tered from unpolarized material, which consists of He, 
15 N, target windows, and foils used to isolate vacuum re- 
gions. The factor Cn — 0.98 corrects for the contribution 
of polarized protons in the 15 N The "pair-symmetric" 
correction C ps and the radiative correction term Arc are 
discussed below. 

The background Nbkg accounts for approximately 85% 
of the detected electrons, and it was determined from an 
analysis of runs using the carbon target. Comparisons 
of the carbon and ammonia scattering rates for W > 
1.4 GeV were used to extract the mass thickness of the 
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ammonia. The known densities and thicknesses of the 
materials in the beam-line, along with a parametrization 
of the ratio of the free proton and neutron cross sections 
to account for the unpaired proton in 15 N, were used to 
calculate Nj,kg from the carbon target data. 

Eq. 0| does not include corrections for detector accep- 
tance or efficiency since, except for possible changes in 
the detector between ammonia and carbon runs, these 
cancel in the ratio. Less than 5% of the data showed 
short-term localized changes in the detector. The af- 
fected regions were removed from the analysis. 

The polarization product P e Pt for each set of runs was 
extracted directly by comparing the measured asymme- 
try for elastic scattering to the known elastic ep asymme- 
try. Two independent procedures were followed. In the 
first method, the background-subtracted asymmetry was 
obtained in several Q 2 bins, checked for consistency, and 
averaged. The second method exploited a limited kine- 
matical range in which both the scattered electron and 
recoil proton could be detected. The strict correlation 
in angle and momentum for elastic ep scattering distin- 
guishes it from inelastic ep scattering and from quasi- 
elastic scattering in nitrogen without subtracting a back- 
ground spectrum. The two methods were consistent with 
each other. A typical value was P e Pt — (36 ± 1)%. 
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FIG. 1: The photon asymmetry A\ + r\A-z vs. the invariant 
hadronic mass W(GeV) for the proton in two Q 2 bins. The 
lower Q 2 data were obtained with a beam energy of 2.6 GeV 
and the higher with 4.3 GeV. The dashed curves show our 
model estimates of r/A2. The error bars denote statistical 
uncertainties and the shaded bands near the bottom of the 
plot indicate the magnitude of the systematic uncertainties 
(la). 

A "pair-symmetric" correction C ps was made for the 
contribution of 7r° — > e + e~7 to the inclusive electron 
sample. It was determined by measuring the e + rate 
with opposite torus current. Since the correction rises 
very sharply at high W (low scattered electron energy), 
we analyzed data only in the kinematic region where the 
pair-symmetric background was less than 10% of the to- 
tal electron rate. 

Radiative corrections were applied to the experimen- 



tal asymmetries using the code RCSLACPOL devel- 
oped at SLAC and based on the approach of Kuchto and 
Shumeiko [2(| for the internal corrections, and by Tsai 
[2?| for external corrections. We used parametrizations 
of the world data on unpolarized and polarized struc- 
ture functions (including our own preliminary asymmetry 
data) and elastic form factors as input for the radiative 
correction code. Details of this model will be given in a 
longer paper. 

The model used for radiative corrections was also used 
to calculate values for A\ , A2, gi , 32 , Fi , Fz, and R over 
the measured region of Q 2 and x. Using the calculated 
value of R, Eq. ^was used to obtain A\ + 77A2 from the 
corrected experimental asymmetry defined by Eq.Q] The 
results for two representative Q 2 bins are shown in Fig.^ 
The model values of 77A2 were generally small, as can be 
seen from the dashed curves in Fig. ^ These calculated 
values were used to extract A\ from the measured A exp . 
Finally, we used the model for F% and A2, along with 
Eq. 21 to extract gi, which is dominated by A\. The 
results for g\ for the proton for five Q 2 bins are shown 
in Fig. El The solid curve shows g\ calculated from the 
data parametrization used in RCSLACPOL. 




FIG. 2: The polarized structure function <?i (x) vs. x for the 
proton for five Q 2 bins. The solid lines show the parametriza- 
tion described in the text. The error bars denote statistical 
uncertainties and the shaded area at the bottom of each plot 
indicates the magnitude of the systematic uncertainties (la). 

We can integrate gi(x, Q 2 ) from the lowest measured 
value of x (with a cut-off ranging from W=2 to 2.6 GeV) 
up to x = 1 at each Q 2 to obtain the contribution of our 
data to the integral T\ (the contribution of elastic scat- 
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tering is excluded). The results are shown as the closed 
circles in Fig. [3] The full value of the integral was esti- 
mated by using the parametrization for gi to estimate the 
contribution from x = 10 -5 up to the threshold of our 
measurements. These extended integrations are shown 
as open circles. The plot is cut off at Q 2 — 1.2 GeV 2 
since at higher values more than half the integral is due 
to the unmeasured region of x. The slope at Q 2 = 
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FIG. 3: The first moment I\ vs. Q 2 for the proton. The 
closed circles are obtained by integrating over our measure- 
ments of gi . The open circles are obtained by using the inclu- 
sive scattering model to extend the integral down to a; = 10~ 5 . 
The error bars indicate statistical uncertainties. The dark 
shaded region at the bottom shows the estimated system- 
atic uncertainties (1 a) for the measured points. The line 
above the shaded region indicates the systematic uncertainty 
for the open circles. Measurements from SLAC |5j are shown 
as shaded squares with the double error bar indicating statis- 
tical uncertainty and (statistical + systematic) uncertainties. 
See the text for an explanation of the curves. 

required by the GDH sum rule is indicated by a straight 
line. Calculations using the light baryon ChPT formula- 
tion of Bernard et al. [I3 are indicated by a shaded band 
at Q 2 < 0.2 GeV 2 with boundaries that reflect uncertain- 
ties in resonance parameters. The pQCD evolution of the 
deep inelastic scattering measurements to O(a^) [l!j is 
shown as a line at high Q 2 . The calculations of Soffer and 
Teryaev j2^| and Burkert and Ioffe [2(| in the intermedi- 
ate region are also shown. A linear fit to the five points 
from Q 2 =0.20 to 0.47 GeV 2 yields a zero crossing for Ti 
at Q 2 = 0.29 ± 0.03 GeV 2 , where the error includes only 
the statistical uncertainty and the estimated uncertainty 
in the DIS contribution added in quadrature. The zero 
crossing indicates the transition to a distance scale where 



non-partonic contributions such as resonance excitations 
are dominant. Our results for Ti(Q 2 ) lie well below the 
predictions from the pQCD evolution from DIS. They 
are in better agreement with the model calculations of 
Ref. [2(J that include s-channel baryon resonance excita- 
tions explicitly. Thus, we think that it is likely that the 
lack of explicit inclusion of the resonance contributions 
in the pQCD evolution gives rise to the discrepancy. 

The estimated systematic uncertainties (1 a) are in- 
dicated by the shaded band at the bottom of each plot 
in the figures. The systematic uncertainty is dominated 
by the parametrizations of A2, f\ and R, which consti- 
tute 75% of the total uncertainty at low Q 2 and 50% 
at high Q 2 . The uncertainty is estimated by using al- 
ternative parametrizations, as well as by setting R and 
A2 to extreme values. The radiative corrections, which 
also incorporate these models, constitute 20% of the to- 
tal systematic uncertainty at low Q 2 and 5% at high Q 2 . 
The remaining uncertainty arises largely from livetime 
calculations and the removal of the nuclear contributions 
from the ammonia spectra. The systematic uncertainty 
of the DIS extrapolation for Fig. [31 was estimated by us- 
ing three different parametrizations for g\ in the low x 
region, Simula |2l| , and the previously mentioned model 
fitted to world data before 1999 and 2000, respectively. 

In summary, our measured asymmetries for inclusive 
scattering of 2.6 and 4.3 GeV polarized electrons on po- 
larized protons have been used to extract the structure 
function 5i(x,Q 2 ) for Q 2 from 0.15 to 1.64 GeV 2 . The 
first moment Ti depends strongly on Q 2 , with its sign 
changing near Q 2 — 0.3 GeV 2 . 
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